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ABSTRACT 

Wc present an analytic formulation for the equilibrium gas density profile of early-type 
galaxies that explicitly includes the contribution of stars in the gravitational potential. 
We build a realistic model for an isolated elliptical galaxy and explore the equilibrium 
gas configurations as a function of multiple parameters. For an assumed central gas 
temperature fesTb = 0.6 keV, we find that neglecting the gravitational effects of stars, 
which can contribute substantially in the innermost regions, leads to an underestimate 
. of the enclosed baryonic gas mass by up to —65% at the effective radius, and by up to 

rj ' ~15% at the NFW scale radius, depending on the stellar baryon fraction. This formula 

. , is therefore important for estimating the baryon fraction in an unbiased fashion. These 

■ new hydrostatic equilibrium solutions, derived for the isothermal and polytropic cases, 

can also be used to generate more realistic initial conditions for simulations of elliptical 
q , galaxies. Moreover, the new formulation is relevant when interpreting X-ray data. We 

5h ' compare our composite isothermal model to the standard /3-model used to fit X-ray 

ryj \ observations of early-type galaxies, to determine the value of the NFW scale radius 

C3 . r s . Assuming a 10% stellar baryon fraction, we find that the exclusion of stars from 

the gravitational potential leads to (i) an underestimate of r s by —80%, and to (ii) 
[ an overestimate of the enclosed dark matter at r s by a factor of ~2, compared to the 

i equivalent /3-model fit results when stars arc not taken into account. For higher stellar 

l/") ' mass fractions, a /3-modcl is unable to accurately reproduce our solution, indicating 

f**". , that when the observed surface brightness profile of an isolated elliptical galaxy is 

found to be well fitted by a /3-model, the stellar mass fraction cannot be much greater 
(N ; than -10%. 
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sj ■ 1 INTRODUCTION 

Several studies, both observational and computational, of groups and clusters of galaxies, and of individual early-type galaxies, 
assume that the hot gas present in these systems is in hydrostatic equilibrium (HE) in the overall gravitational potential. 
Even though the assumption of HE has been shown to not always be robust (e.g. Bertin et al. 1993; Ciotti & Pellegrini 2004), 
many authors assume an isothermal /3-model (Cavaliere & Fusco-Femiano 1976, 1978) to describe the gas profile, but several 
solve the HE equation, assuming that the gravitational potential is due to dark matter (DM) only, usually described by an 
NFW (Navarro, Frenk & White 1996) profile (e.g. Makino, Sasaki & Suto 1998; Suto, Sasaki & Makino 1998), or by more 
general profiles (e.g. Ciotti & Pellegrini 2008). While this approximation seems to hold for groups and clusters of galaxies, 
it is well known that stars contribute a significant fraction of the total and baryonic matter of elliptical galaxies, and even 
become the dominant part at the very centre of these systems (e.g. Ferreras, Saha & Williams 2005). The purpose of this 
paper is to extend previous calculations that only considered the potential due to DM, and explicitly include the contribution 
of the stellar component as well. 

Providing an analytic fitting formula for calculating the response of the gas to the presence of the gravitational potential 
induced by both DM and stars is useful while interpreting X-ray data of early-type galaxies and for setting up initial conditions 
for simulations of elliptical galaxies. The derivation of an equilibrium gas profile neglecting the effects of stars leads to an 
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incorrect estimate of the enclosed gas mass, particularly in the central regions of the galaxy, and therefore has implications 
for the evaluation of the baryon fraction as a function of radius, and for the calculation of the rate of accretion onto a central 
supermassive black hole. Moreover, the known relations (e.g. Makino et al. 1998) between observed properties of the hot gas, 
described via a /3-model, and the parameters of the underlying DM gravitational potential, are significantly affected by the 
inclusion of the stellar component. In addition, an analytic formula that includes the effects of stars can be important to 
provide a more realistic formulation of the initial equilibrium conditions for simulations of galaxy formation and evolution 
(e.g. Di Matteo, Springel & Hernquist 2005; Springel, Di Matteo & Hernquist 2005; Hopkins et al. 2006). 

It has long been observed, both at low redshift (e.g. Graham et al. 1996; Kormendy et al. 2009) and high redshift (e.g. 
Mancini et al. 2010; van Dokkum et al. 2010), that the surface brightness of many elliptical galaxies can be well described 
by the empirically derived de Vaucouleurs (1948) profile or, more generally, by the Sersic (1963, 1968) profile. Unfortunately, 
the de-projected stellar mass density profile cannot be calculated analytically (but see Young 1976 for a numerical study). 
There is, however, a family of stellar mass density profiles (Dehnen 1993; see also Tremaine et al. 1994) that also describes 
observed properties well. Not only is the Dehnen family analytically simple, but also, for a defined range of parameters, it 
quite successfully resembles the de Vaucouleurs profile in both surface density and distribution function. Furthermore, it 
includes, as special cases, the widely used Jaffe (1983) and Hernquist (1990) profiles. It is this family of density profiles that 
we investigate further in this paper, and use to model the stellar component. Throughout the paper, we use the following 
values for cosmological parameters when needed: fio = 0.27, fit, = 0.045, J7a = 0.73, Ho = 100 h km s _1 Mpc -1 = 70/i7o km 
s" 1 Mpc" 1 , h = 0.7 (i.e. h 70 = 1.0). 

The outline of the paper is as follows. In Section [2] we formulate the problem and derive the new formalism for HE gas 
profiles in the presence of DM and stars. In Section [3] we analyse the HE solutions for a realistic model of an isolated elliptical 
galaxy, compare the results for a wide range of parameters, and present new relations between /3-model and DM parameters. 
We conclude in Section [4] with a discussion of our results and their implications for interpreting X-ray observations, and for 
simulations of elliptical galaxies. 



2 FORMULATION 

We model the composite system of gas in equilibrium with stars and a DM halo. To start with, we consider the case of an 
ideal gas in a generic total gravitational potential <f) with a finite central value <j>o, and solve the HE equation, VP = —pV(f>, 
where P and p refer to the pressure and mass density of the gas. Note that the general equation simplifies in the case of 
spherical symmetry to the familiar form 



dP 
dr 



GM T (r) 



(1) 



where we have made use of Poisson's equation, V 2 </> = 4rcGpT, and Mt(t) and pr(r) are the total enclosed mass and total 
mass density, respectively. If the gas is barotropic, i.e. p — p(P), we can further write 
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where the free parameter Po — P(0) was chosen as the boundary condition. 

For simplicity, we only consider the two common barotropic cases: the isothermal and the polytropic case. In the isothermal 
ideal gas case, pressure and density are related by P(r) = [kBTo/(pm p )]p(r) — (Po/po )/?(»")> where fcs is the Boltzmann 
constant, p is the gas mean molecular weight, m p is the proton mass, and To and po are the central gas temperature and 
density, respectively. The solution to the HE equation in this case can be reduced to the familiar integral: 
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where A = —<j>opo/Po- 

For the polytropic ideal gas case, pressure and density are related by P(r) = Kop(r) r , where Ko 
fcsTb/(A tm pPo _1 ) ano - T 7^ 1 is the polytropic index. The solution to the HE equation in this case is: 
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(4) 



It is worth pointing out that the polytropic solution is only physically meaningful if the term in square brackets of the 
r.h.s of equation Q is non-negative. Moreover, the isothermal quantities (and, in most cases, the polytropic quantities) are 
non-zero at infinity: with our choice of the boundary condition, regardless of the slope of the gravitational potential at large 
radii, only a potential which is infinite at the origin would be able to cause the gas density at infinity to be zero. For this 
reason, results at large radii should be treated with caution. See Bulbul et al. (2009) for an alternate approach in the polytropic 
case, where they impose T(+oo) = as the boundary condition. This is equivalent to setting (F— 1)A/F = 1 in equation Q. 
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We further assume that gas does not contribute to the total gravitational potential (i.e. we neglect gas self-gravity), and 
therefore solve the simplified HE equation instead, where now <j> is the gravitational potential due to all mass excluding gas, 
and Mt is the total mass excluding gas. This assumption is supported, especially for the inner regions, by previous studies 
(e.g. Suto et al. 1998), and is justified by the fact that gas is globally negligible, i.e. the enclosed gas mass is one or more 
orders of magnitude smaller than the total enclosed mass for all radii up to the virial radius. Stars, on the other hand, while 
having a virial enclosed mass comparable to that of the gas, are not negligible in the inner parts of the galaxy: in fact, they 
represent the dominant matter fraction at least up to the effective radius (e.g. Ferreras et al. 2005). We initially consider a 
system with gas and DM only, and assume that the DM density can be described by an NFW profile, 



pNFw{r) = . . . j— — 7— rj, (5 
(r/r s )(l + r/r s y 

where r 3 is a scale radius, 8 C is a characteristic (dimensionless) density, and p c is the critical density of the Universe. The 
corresponding gravitational potential and enclosed mass are 



4>NFw(r) = 4>nfwo — i}—tl/lfl M NF w(r) = 47r<5 c p c rJ 
r/r s 



In | 1 + — 



(6) 



respectively, where (Pnfwo = (Pnfw(0) = —4TvGS c p c r 3 , and we have assumed 4>nfw(+oo) — 0. Substituting these into 
equations © and @, we have, as in Makino et al. (1998) for the isothermal case, 



p(r) = p exp 



-Anfw ( 1 



ln(l + r/r 3 ) 
r/r 3 



and, as in Suto et al. (1998) for the polytropic case, 
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(8) 



where A NF w = -<PnfwoPo/Po = 4TTGS c p c r 3 pm p /(kBT ). 

Now we move on to the case of an NFW DM halo with the stellar component modelled using a Dehnen profile, neglecting 
the effects of adiabatic contraction. The Dehnen (1993; see also Tremaine et al. 1994) family of stellar mass density profiles 
is given by 



(3 - t/QAf, r. 

PD{r)= 4tt r*(r + r.) 4 -*' (9) 

where < ip < 3, M* is the total stellar mass (i.e. integrating from the origin to infinity), and r* is a scale radius. The Dehnen 
profile is useful as its projected distribution accurately matches the de Vaucouleurs (1948) surface density profile which is an 
excellent fit to observations. This is especially true for values of ip ranging between 1 (which corresponds to the Hernquist 
profile) and 2 (the Jaffe profile), with ip = 3/2 being the value for which the match is closest. The Dehnen density profile 
corresponds to a gravitational potential 



4>n(r) = 



GM, 



2-ip 



r + r* 



if V + 2, 
if V = 2, 



(10) 



where we have assumed </>.d(+oo) = 0. Since equations ((3J and (j4| require a finite central gravitational potential, we only 
consider Dehnen profiles with < tp < 2, for which the gravitational potential and enclosed mass are 



M D (r) = M^-^-J , (11) 

respectively, where 4>do = 0n(O) = -GM»/[r,(2 — tp)]. 

We consider a system that consists of an NFW DM halo and stars described by a Dehnen profile, for which the total 
gravitational potential and total enclosed mass are given by explicitly adding the contributions of the stars and DM given in 
equations (JTTJ) and l[6]l. respectively. Using equations © and (Q, we obtain 
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for the isothermal case, and 
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(13) 



for the polytropic case, where Ad = —<f>DoPo/Po = GM*/j,m p /[r t kBTo(2 — </))]. Notice that equations (|12H13[) reduce to 
equations (|7I8[) when M* = 0. 



3 ANALYSIS 

3.1 Exploring parameters and realistic matter distributions 

In this section, for a system with an NFW DM halo and stars described by a Dehnen profile, we compare equilibrium gas 
density profiles as a function of f s tar, the stellar mass fraction at the virial radius, t/>, the Dehnen parameter, and To, the 
central gas temperature. We also consider the case when stars are included in the total mass budget, but not in the total 
gravitational budget. 

We consider a system at redshift z, which consists of an isolated elliptical galaxy at the centre of an isolated, collapsed 
DM halo, wherein the density is given by 

PNFw(r) = ; -7— — rj, (14) 

47r j(c) r(r + r s ) 

where MjvFw( r w) is the DM mass enclosed within the virial radius r vir , and f(c v i r ) = ln(l + c vir ) — c vir /(l + c vir ) is a 
function of the concentration parameter cw = r v i r /r s , which depends on redshift and (total) virial mass M v i r as 

_ cp / M vir \ a 

Cmr ~ 1 + z \1.3xW^h-^Me) ' 1 j 

where Co and a are constants that can be inferred from simulations (e.g. Bullock et al. 2001; Hennawi et al. 2007) and 
observations (e.g. Comerford & Natarajan 2007; Mandelbaum, Seljak & Hirata 2008), and M V i r is in units of Mg. 

The virial radius is, by definition, the radius of a spherical region within which the mean (total) mass density is A V i r 
times the mean critical densitjQ: 

/ W ■ \ 1//3 

(!6) 



' vir — I . . 

where the virial overdensity A vir (z) ~ 187r 2 + 82(Q m (z) — 1) — 39(f2, n (z) — l) 2 (Bryan & Norman 19980, the matter fraction 
ttm(z) = fio(l + z) 3 /[n (l + z) z + tt A ], and the mean critical density p c (z) = 3J?oPo(l + z) 3 + Qa]/(87tG). 

The NFW density profile function can be also re-written as in equation ([5]), with 5 C — /cA/A 1 ji r c 3 ir /(3/(c l ,i r )), and 
fr>M = Mnfw (rvir) /M V ir = 1 — fb = 1 ~ bf,Slb/fio, where /dm and are the DM and baryon fraction, respectively, within 
the virial radius, bb is the baryon fraction relative to the universal value, and we have assumed that the universal baryon 
fraction, Qj,/f2o, is independent of redshift. 

The elliptical galaxy stellar density profile is given by equation The total stellar mass M* is calculated by setting 
r = r V i r in equation (|ll|l and by imposing AlD(r V ir) = fstarM v i r . The virial stellar mass fraction can be written, following 
e.g. Mamon & Lokas (2005b), as 

J star — „ r , . — — , (II ) 

I B-L*,B(r v ir) OTIS 

where T^b is the stellar mass-to-light ratio, T B is the universal mass-to-light ratio, Ts is the galactic mass-to-light ratio 
(all in M0/L0), and &y is the (dimensionless) mass-to- light ratio bias, all in the B-band. The scale radius r* is related to the 
effective radius r e through a numerical fit (Dehnen 1993), 



2 i/(3-V0 



r e 0.7549 - 0.00439-0 + 0.00322^ 2 - O.OO1820 3 



(18) 



1 We use p c {z); other authors, e.g. NFW, use p c (0). 

2 With the above definitions, A„ ir (0) ~ 97, A„ ir (0.5) ~ 133, and A„ iT .(+oo) = 18tt 2 . Note that other authors (e.g. Bullock et al. 2001) 
use equivalent but different definitions, M v i r = 4irr 3 ir A v i r Q m (z)p c /3 and A v i r (z) ~ [18-7T 2 + 82(f7 m (z) — 1) — 39(f2 m (z) — l) 2 ]/Q m (z), 
so that A v i r (0) ~ 359. 

3 Notice that by imposing the same stellar virial mass (i.e. by fixing fstar) and by varying the Dehnen parameter, the total stellar mass 
M, will necessarily vary, albeit only very slightly. 
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Cosmologies! parameters 


o„ n 97 o . n 71. o n n/t ^ v, n 7 h 1 

lio — u.z / , i — u./o, ii^ — u.u^to, n — u. / , Ofj — i 


Stellar parameters 


Of — U.ZOD'il, 1 b — oyU IVIq / 1j0 , 1 # B — 0.0 IVIq / Ijq , 1/J — 1 


Gas parameters 


1 — U.D Ke V,/J — U.DZ, 1 — / O 


Virial mass and redshift 


M„ 8r = 10 12 M , z = 0.5 


Concentration parameters 


c = 9, q = -0.13 


Gas central density 


po = 5.80 x 10" 27 g cm" 3 


Gas central pressure 


Po = 5.37 x 10~ 12 g cm" 1 s" 2 


Matter fraction, critical density, and virial over- 
density at redshift z 


ft z = 0.56, p c = 1.51 x 10" 29 g cm" 3 , A„„. = 133 


Virial radius and NFW parameters 


r v „ = 200 kpc, r s = 22.8 kpc, c mT = 8.77, S c = 1.81 x 10 4 


Central gravitational potentials 


(po-NFW = -1.14 x 10 15 cm 2 s~ 2 , cfo-D = -2.35 x 10 15 cm 2 s~ 2 


Shape parameters 


Anfw — 1.23, — 2.53 


Stellar scale radii 


r e = 2.19 kpc, r* = 1.21 kpc 



Table 1. Relevant parameters used in our fiducial model, for an isothermal gas density profile. Given quantities are in grey cells; 
calculated quantities are in white cells. 



whereas the effective radius itself can be derived from empirical fits to observation^ (Mamon & Lokas 2005a) as log(hror e ) = 
0.34 + 0.54 log Lio + 0.25(log Lio) 2 , where Lio = h? L B /{W 10 L Q ), L B = M vir /{brT B ), and L B and M vir are in units of L© 
and Mq, respectively. 

Finally, we can calculate the central gas density po by imposing 

/ p(r)Anr 2 dr = f gas M v i r = (1 - fstar - fDu)M vir = ( b b ^- - t£ j M vir , (19) 
Jo V iZ o OtTb/ 

where we have assumed that all matter is comprised of DM, gas, and stars only (i.e. we have neglected other forms of matter, 
such as dust, and the presence of a supermassive central black hole), and that p(r) is given by equation (|12[) or (|13[) when 
stars are included in the gravitational potential, and by equation ([7]) or ([8J when stars are excluded. Note, however, that the 
gas virial mass is the same with and without the stellar component in the gravitational potential, since stars are included in 
the total mass budget in both cases. 

We note that a very small number of parameters are needed to compute the equilibrium profiles for the gas. These are, 
for the DM halo system, the redshift z, the (total) virial mass M V i r of the system, the fitting parameters Co and a of the 
empirical relation between concentration and total virial mass, and the baryon fraction bb relative to the universal value. For 
the stars, the parameters are the galactic B-band mass-to-light ratio T B (or, alternatively, the universal B-band mass-to-light 
ratio Ts and the B-band mass-to-light ratio bias br), the stellar B-band mass-to-light ratio T*,b, and the Dehnen parameter 
ip. Finally, for the gas, the parameters are the mean molecular weight p,, the polytropic index F, and the central temperature 
To (or, alternatively, the ratio between central density and pressure pa /Pa). 

In our fiducial model, we use the Mamon & Lokas (2005b) mass-to-light ratio values and choose &r = 0.25641, T B = 390 
Mq/L© (so that T_b = 100 Mq/L© precisely), and T*,s = 6.5 M0/L0. As for the halo parameters, we choose a system 
with universal baryon fraction (bb = 1) at redshift 2 = 0.5, a (total) virial mass M V i r = 10 12 Mq, and the c V i r -M V i r relation 
parameters co = 9 and a = —0.13 (Bullock et al. 2001). Finally, we fix the central temperature of the gas to be k B To = 0.6 
keV (Diehl & Statler 2008), the mean molecular weight of the gas to be p = 0.62, the polytropic index to be F = 5/3, and 
the Dehnen parameter to be ip = 1 (Hernquist 1990). All these quantities are listed in Table [U in addition to the parameters 
calculated throughout this section. 

The choice of the central temperature Tq is particularly important, since it greatly influences the values of Anfw and 
Ad. If we use the (cluster and group) mass-temperature (M-T) relation (e.g. Sanderson et al. 2003) or, rather, extrapolate 
the M-T relation down to the mass ranges typical of early-type galaxies, we obtain e.g. k B To — 0.2 keV for a system with 
M v i r = 10 12 Mq. However, the extrapolation of the M-T relation is not without risk, as there is some evidence for a gradual 
steepening of the relation, with decreasing mass (Sanderson et al. 2003). Using individual observations of isolated elliptical 
galaxies, the central temperatures measured are usually slightly higher, ranging from k B To ~ 0.5 keV (e.g. Memola et al. 
2009) to k B T ~ 1 keV (e.g. O'Sullivan & Ponman 2004; O'Sullivan, Sanderson & Ponman 2007). Finally, Diehl & Statler 
(2008) give a convincing argument for a central temperature k B To ~ 0.6 keV, using observations of 36 normal early-type 
galaxies. For this reason, here we study cases with three different central temperatures: k B To = 0.2, 0.6, 1 keV. 

In Fig. [T]we plot the density (top panels) and enclosed mass (middle panels) profiles of DM, stars, and isothermal gas, for 
the fiducial model (with k B To = 0.6 keV; central panels) and for the lower and higher temperature models (with k B To = 0.2 
keV, left panels; and k B To = 1 keV, right panels), with and without the stellar component in the gravitational potential. We 



4 Notice that r e does not depend on the stellar fraction. 
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k B T = 0.2 keV k B T = 0.6 keV k B T = 1.0 keV 




-3-2-10 1 2 -3-2-10 1 2 3 -2-10 1 2 3 

'og, ( r A s ) 

Figure 1. Mass density, enclosed mass, and fractional gas mass change as a function of radius, for an elliptical galaxy with DM (NFW 
profile), stars (Hernquist profile), and isothermal gas. For the fiducial model (central panels), whose parameters are listed in Table[l] and 
for two models with lower (left panels) and higher (right panels) gas temperature, we plot in the top (middle) panels, the mass density 
(enclosed mass) profiles of DM (red, dotted curve), stars (red, dashed curve), and isothermal gas (blue, dot-dashed curve) in HE with 
DM and stars. The black, solid curve represents the isothermal gas in HE with DM only. In the bottom panels, we plot (solid curve) the 
fractional change of enclosed gas mass between the case when stars are included in the total gravitational potential, and the case when 
they are not; we also plot (dashed curve) the ratio between the change in gravitational energy of a gas particle, due to the addition of 
the stellar component, and its thermal energy. The four vertical lines denote the stellar scale radius r», the effective radius r e , the NFW 
scale radius r s (solid line), and the virial radius r„; r , from left to right, respectively. The inclusion of the stellar component dramatically 
alters the gas profile at the very centre of the early-type galaxy, with the gas being more centrally concentrated, due to a deeper and 
steeper gravitational potential well. The effect of gas concentration is more pronounced for lower values of the gas temperature. 

focused on the isothermal solution because the radial range for which the polytropic solution is non-negative is limiteojfl. The 
inclusion of the stellar component dramatically alters the equilibrium profiles for the gas at the very centre of an early-type 
galaxy, with the gas being more centrally concentrated, due to a deeper and steeper gravitational potential well. The presence 
of the stellar component is thus explicitly taken into account and is manifested in the modification to the gas profile in the 
central region, in contrast to the case when stars are not included in the total gravitational potential. In the bottom panels 
of Fig. [TJ we plot the fractional change of enclosed gas mass AM/M as a function of radius: for the fiducial model, in the 
central panel, the gas concentration ratio between the new (i.e. including stars in the gravitational potential) and the old (i.e. 
excluding stars) enclosed gas mass is [1 + AM/M] = f CO nc ~ 4 at the stellar scale radius r», ~ 3 at the effective radius r e , 
and ~ 1.2 at the NFW scale radiufl r s . We note that the effect of gas concentration increases with decreasing temperature, 
because the ratio between the change in gravitational energy of a gas particle, due to the addition of the stellar component, 
and its thermal energy, decreases with increasing temperature, as is shown in the bottom panels of Fig. [1] We also point 
out that our gas (density and enclosed mass) profile is very different from that found in Mamon & Lokas (2005b). They 



5 With the chosen parameters, the threshold temperature value above which the polytropic solution is physically meaningful at all radii 
is k B T ~ 1 keV. 

6 The fit for the effective radius can be considered uncertain to a factor of 2 (Mamon & Lokas 2005a) : for this reason, we multiplied and 
divided r e by 2, fixing all other parameters, and found that the effect of gas concentration increases as the effective radius decreases, as 
expected. 



© 0000 RAS, MNRAS 000, 000-000 



Hydrostatic equilibrium profiles in elliptical galaxies 7 
T., B =6.5 M /L Q T., B =4.0 M /L T., B =1.5 M /L 




-3-2-10 1 2 -3-2-10 1 2 3 -2-10 1 2 3 

'og, ( r /0 



Figure 2. Mass density, enclosed mass, and fractional gas mass change as a function of radius, for an elliptical galaxy with DM (NFW 
profile), stars (Dehnen profiles), and isothermal gas. For the fiducial model (with T„ g = 6.5 Mq/Lq, fstar/fb — 0.39; left panels), 
whose parameters are listed in Table [Tl and for two models with lower stellar baryon fraction (T» g = 4.0 Mq/Lq, f s tar/fb — 0.24, 
central panels; and T» g = 1.5 Mq/Lq, fstar/fb — 0.09, right panels), we plot in the top (middle) panels, the mass density (enclosed 
mass) profiles of DM (red, dotted curve), stars, and isothermal gas in HE with DM and stars. The stellar (gas) profiles are the three 
upper, red (lower, blue) bundled curves, for three values of the Dehnen parameter: tp = 0.5 (dashed curve), -0 = 1 (solid curve), and 
i/j = 1.5 (dot-dashed curve). The black, long-dashed curve represents the gas profile when the stellar component of the total gravitational 
potential is neglected. In the bottom panels, we plot the fractional change of enclosed gas mass between the case when stars are included 
in the total gravitational potential, and the case when they are not, for the same three values of the Dehnen parameter. The three 
vertical lines in each panel denote the effective radius r e , the NFW scale radius r e (solid line), and the virial radius r v i r , from left to 
right, respectively. The effect of gas concentration is more pronounced for higher values of the stellar mass fraction and for higher values 
of the Dehnen parameter. 



assume a fixed /3-model to describe the isothermal gas in the galaxy, and do not solve the HE equation directly. Note that, 
by construction, the enclosed gas mass at the virial radius is the same. 

It is worth pointing out that the value of fstar of the fiducial model corresponds to a value of stellar baryon fraction 
fstar/fb — 0.39, which is higher than that inferred from recent gravitational lensing studies, for which fstar/fb — 0.1 (e.g. 
Hoekstra et al. 2005; Mandelbaum et al. 2006; Heymans et al. 2006; Lagattuta et al. 2009; see also Fukugita, Hogan & Peebles 
1998). The discrepancy is likely due to different values of mass-to- light ratios and of mass-to-light ratio bias. Lagattuta et 
al. (2009), for example, use a V-band virial mass-to-light ratio Ty = 210 Mq/Lq and a V-band stellar mass-to-light ratio 
Y*,y = 2.8 Mq/Lq. Using these values, we would obtain fstar/fb — 0.1. For this reason, we compare the gas radial profiles 
for three different values of Y,,s = 1.5, 4.0, 6.5 Mq/Lq, which correspond to stellar mass fractions fstar — 0.015, 0.04, 0.065, 
and to stellar baryon fractions fstar/fb — 0.09,0.24,0.39, respectively. For these values of Y*,b, we also vary the Dehnen 
parameter to explore dependencies on stellar profiles, choosing values of ?/> = 0.5, 1 (Hernquist 1990), and 1.5 (the "best" 
value in Dehnen 1993). 

Fig. [2] shows the density (top panels) and enclosed mass (middle panels) profiles for the three chosen values of T,^ and 
of tp, in the same way we plotted Fig.[TJ for a system with kgTo = 0.6 keV. The effect of gas concentration is obviously more 
pronounced for higher values of the stellar mass fraction and for higher values of the Dehnen parameter, as it is shown in the 
bottom panels. However, the variation of the Dehnen parameter does not cause significant changes in the enclosed gas mass, 
except in the very inner-most regions of the galaxy (r < r e /10): for example, for the T^s = 6.5 Mq/Lq case, the enclosed 
gas mass at the effective radius varies by less than 1% when changing the Dehnen parameter from ip = 1 to ip — 0.5, or to 
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Figure 3. For an elliptical galaxy with DM (NFW profile), stars (Hernquist profile), and isothermal gas, we plot the equilibrium gas 
density profiles excluding (top panel) and including (bottom panel) the effects of stars on the total gravitational potential, for a B-band 
stellar mass-to-light ratio T, g = 1.5 Mq/Lq (fstar/ fb — 0.09), for two different gas temperatures: kgTo = 0.6 keV (black, solid lines) 
and kgTo = 1 keV (black, dashed lines). The dotted, red curves are the best-fitting ,3-models in the 0.01 r s < r < 10 r s range, denoted 
by the two dotted vertical lines. The inclusion of stars significantly affects the relation between the /3-model scale radius and the NFW 
scale radius. For higher stellar mass fractions (not shown here), a /3-model is unable to accurately reproduce our solution. When the 
observed surface brightness profile of an isolated elliptical galaxy is found to be well fitted by a /3-model, the stellar mass fraction cannot 
be much greater than ~10%. 



tp = 1.5. On the other hand, the change in stellar fraction is much more evident: for example, for the Hernquist case (ip = 1), 
the ratio between the enclosed gas masses at the effective radius r e varies from ~3 in the T,,s = 6.5 M©/L<0 case, to ~2 in 
the Y*,b = 4.0 M0/L0 case, to ~1.3 in the T*,b = 1-5 Mq/Lq case. 



3.2 Comparison with the /3-model 

By comparing the theoretically derived equilibrium solution for an isothermal gas to the observationally inspired /3-model 
profile (Cavaliere & Fusco-Femiano 1976, 1978), it is possible to derive relations between observed properties of the hot gas 
and the parameters of the underlying DM gravitational potential of a given system. Makino et al. pursue this (1998; see also 
Suto et al. 1998 and Wu & Xue 2000 for other DM profiles) for the case of clusters of galaxies with a DM halo described by 
an NFW profile. They fit their HE analytic solution for an isothermal gas (our equation [7]) to a generic /3-model profile given 
by 

( r \ _ PQ3 = PojAp + A 1 A NF w) 

PPV> [1 + (r/r c ) 2 p/ 2 [1 + [r/{A 2 r a )ff A ^NFw^ ' yv > 

where pop, r c , and /3 are the central density, the scale radius, and the characteristic slope parameter of the ,3-model, respectively, 
and Ai (for i = 0-3) are fitting parameters, in the 0.01 r s < r < 10r s range, obtaining r c ~ 0.22 r s and /3 ~ 0.067 Anfw, for a 
wide range of values of Anfw- We performed a non- linear least squares fitting with the IDL routine MPFIT (Markwardt 2009), 
using our own values of A^pw, and obtained similar results in the case when stars are not included in the total gravitational 
potential, regardless of the value of the stellar mass fraction and of the gas temperature: r c ~ 0.239 r a and /3 ~ 0.067 Aatfw 
in the 0.01 r s < r < 10 r s range. Note that, even though the /3-Anfw relation is independent of temperature, the value of 
the characteristic slope parameter /3 does depend on temperature through Anfw- 

Now we go a step further and compare the /3-model profile to the isothermal HE gas profile in the case when stars are 
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included in the gravitational potential, given by equation (|12[) . Since we found that there is no strong dependence on the 
Dehnen parameter, we focus on the Hernquist case (if) = 1) for simplicity. Also, we choose the stellar B-band mass-to-light ratio 
to be T»,b = 1.5 Mq/L©, as the inferred stellar baryon fraction fstar/fb — 0.09 is closer to what found from gravitational 
lensing observations (e.g. Hoekstra et al. 2005), and because higher stellar mass fraction systems do not fit the /3-model 
profiles as well. The fact that it is indeed possible to fit X-ray observations of early-type galaxies with /3-model profiles (e.g. 
O'Sullivan, Ponman & Collins 2003), is an additional indication that the stellar baryon fraction of such galaxies cannot be 
much higher than — 10%. 

In Fig. [3] we show the equilibrium gas density profiles for an isothermal gas in a system with NFW DM and Hernquist 
stars, for two different temperatures, including and excluding the effects of stars on the total gravitational potential. All HE 
gas profiles are fitted to /3-model profiles in the 0.01 r s < r < 10 r 3 range. We chose this particular range to be able to make 
consistent comparisons with the results of Makino et al. (1998), and also because it is typically difficult to obtain precise X-ray 
measurements at larger radii (due to the low surface brightness of the gas) and at smaller radii (due to resolution limitations). 
We note that, whereas the relation between /3 and Anfw does not change significantly (/3 ~ 0.067 Anfw excluding stars, 
and j3 ~ 0.069 Anfw including stars), the r c -r s relation between the two scale radii varies substantially, from r c ~ 0.239 r s , 
when stars are not included in the potential, to r c ~ 0.053 r 3 , when stars are included. We note that the particular constant of 
proportionality between r c and r s depends on the assumed stellar mass fraction. Using the r c -r s relation can therefore lead to 
an underestimate of the NFW scale radius, when stars are not taken into account. Defining f rs = r s - n ew /r s ~oid as the ratio 
between the /3-model inferred value of r a when stars are included in the gravitational potential (the new value), and when they 
are not (the old value), we have f r „ ~ 5. Assuming the virial mass M„, r and radius r v i r of the system are known, we then have 
the ratio of concentration parameters c V i r -new /cvir-oid — l//r„ ~ 0.2. Therefore, the enclosed DM mass at r a - new would be 

Overestimated by a factor Of [ln( 1~\~ fr s Cyir — new ) fr s Cvir — new j fr s C-vir — new )] / [ln(l - ( - Cujr — new ) Cy ir — new j ( l~t~Ci; jr — new )] r * J 2. 

4 CONCLUSIONS 

We solve the HE equation analytically for a system with DM described by an NFW profile, and stars described by a Dehnen 
profile, for the two common barotropic cases (the isothermal case and the polytropic case). We thus extend previous calculations 
which only considered the potential due to DM. In doing so, we do assume spherical symmetry, and we neglect the effects 
of gas self-gravity and adiabatic contraction. We compare the newly obtained equilibrium solutions as a function of stellar 
mass fraction, central gas temperature, and the Dehnen parameter, to the case when the stellar component is included in the 
total mass budget but not in the total gravitational potential. We determine that neglecting the effects of stars leads to an 
underestimate of the enclosed gas mass in the inner regions of the galaxy. We then present new relations between /3-model 
and DM parameters, and find that the inclusion of the stellar component is necessary to avoid an underestimate of the NFW 
scale radius. In particular: 

• For a system with NFW DM and Hernquist stars at redshift z = 0.5, with virial mass M v i r = 10 12 Mq, and central 
temperature ksTo = 0.6 keV, we find in the isothermal case that the gas concentration ratio f CO nc between the new (i.e. 
including stars in the gravitational potential) and the old (i.e. excluding stars) enclosed gas mass, can be as high as ~3 at 
the effective radius, and as high as — 1.2 at the NFW scale radius, for a stellar baryon fraction fstar/fb — 0.39. 

• The effect of gas concentration decreases with decreasing stellar baryon fraction as expected, with f concise) varying from 
—3 in the fstar/fb — 0.39 case, to —2 in the fstar/fb — 0.24 case, to — 1.3 in the fstar/fb — 0.09 case. 

• On the other hand, the effect of gas concentration is virtually independent of the Dehnen parameter, except in the very 
inner-most regions of the galaxy (r < r e /10), with the value of f CO nc(re) changing by — 1% when varying the Dehnen parameter 
in the 0.5 < ip < 1.5 range. 

• We further determine that the effect of gas concentration decreases with increasing temperature, as expected. 

• We find similar results for the polytropic case at higher temperatures, where the solution is physically meaningful. In fact, 
for a fixed central gas temperature, the effect of gas concentration is more pronounced than in the isothermal case, because 
the polytropic gas has a lower mean temperature. 

• We derive new relations between observed properties of the hot gas, described via a /3-model, and the parameters of 
the underlying DM gravitational potential of early-type galaxies, which avoid an underestimate of the NFW scale radius by 
—80%, and an overestimate of the enclosed DM mass at r s by a factor of —2, assuming a 10% stellar baryon fraction. 

• Models with higher stellar mass fraction do not fit /3-models well, further indicating that the stellar baryon fraction of 
isolated elliptical galaxies whose observed surface brightness profile is well fitted by a /3-model cannot be much higher than 
-10%. 

These new analytic solutions will be useful for the interpretation of X-ray data of elliptical galaxies and in particular, 
will help derive unbiased estimates of the baryon mass fraction as a function of radius. Moreover, the formalism in equations 
(|12H13p provides a new prescription for setting up more realistic initial conditions for simulations of galaxy formation (e.g. Di 
Matteo et al. 2005; Springel et al. 2005; Hopkins et al. 2006). 
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